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Molecular fluids show rich and complicated dynamics close to the glass transition. Some of these 
observations are related to the fact that translational and orientational degrees of freedom couple 
in nontrivial ways. A model system which can serve as a paradigm to understand these couplings 
is a system of hard ellipsoids of revolution. To test this we compare at the ideal glass transition 
the static molecular correlators of a linear A-B Lennard-Jones molecule obtained from a molecular 
dynamics simulation with a selected fluid of hard ellipsoids for which the static correlators have 
been obtained using Percus-Yevick theory. We also demonstrate that the critical non-ergodicity 
parameters obtained from molecular mode coupling theory for both systems show a remarkable 
similarity at the glass transition, provided the aspect ratio is chosen properly. Therefore we conclude 
that a system of hard ellipsoids can indeed be used to understand part of the essential behaviour 
of such a simple molecular system like the A-B Lennard-Jones molecules in the vicinity of the ideal 
glass transition. 
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I. INTRODUCTION 

When a fluid is driven towards a glass transition ei- 
ther by quenching down the temperature (super-cooling) 
or by increasing the density a critical slowing down of the 
dynamics occurs without a diverging length scale in the 
system. 

However there are only a few systems where the full dy- 
namics of the glass transition is understood in detail and 
where it is possible to make predictions for the criti- 
cal density and temperature for the ideal glass transi- 
tion from a microscopic theory. Systems where such a 
transition is well understood are colloidal suspensions |Q . 
Their liquid state is well described by a hard sphere sys- 
tem and the static structure factor calculated by Percus- 
Yevick approximation has been used as an input for the 
mode coupling theory (MCT) Q to calculate the dynam- 
ics near the ideal glass transition ^Jjl. Within idealized 
MCT this transition takes place at a critical tempera- 
ture T c or a critical density p c at which a breaking of 
ergodicity occurs. For a simple glass former like a col- 
loidal system the mode coupling equations describe the 
dynamics in great detail, explain why the transition is 
dominated by the cage effect and their solutions are in 
good agreement with experimental data . 
Most glass formers do not have such a simple structure 
like hard spheres. A real glass former is usually a molec- 
ular system and shows a wide variety of additional struc- 
ture in the dynamics (multiple microscopic peaks, fast 
and slow /3 peaks, a, a' peaks) where a detailed under- 
standing is lacking so far. A large part of this additional 
structure is caused by the coupling of translational and 
orientational degrees of freedom. There are different pos- 
sibilities to get a deeper insight into such couplings. 
The hrst possibility to study the relaxation in the super- 
cooled regime is to model a specific glass former as ac- 
curately as possible using model potentials, and then to 



perform computer simulations. This has been done, e.g. 
for water for silicate glass |t| or for a linear A-B- 

molecule [pf. 

A second possibility to get insight into the statics as well 
as in the dynamics of a glass former is to study model 
systems which are as simple as possible but keep the es- 
sential mechanism - in our case the interplay between 
translational and rotational degrees of freedom. Such a 
model system which can be used is a system of hard el- 
lipsoids of revolution. 

The aim of this paper is to answer the question how far 
a system of hard ellipsoids can help to understand basic 
mechanisms present in molecular glass formers. 
We therefore compare in detail the static structure of 
the translational and orientational degrees of freedom 
obtained from a molecular dynamics simulation for an 
A-B-Lennard Jones molecule with results for hard ellip- 
soids of revolution. Using the static correlators as an 
input into the mode coupling equations we also compare 
the results for the critical non-ergodicity parameters with 
each other. 



II. MODE-COUPLING EQUATIONS 

As a first step towards finally calculating the full dy- 
namics we want to test for both systems the occurrence 
of ergodicity breaking. Therefore the MCT- equations 
are solved in the limit t — ► oo The non-ergodicity 
parameters are given by (F(q,m)) u , = Fu'(q,m) = 
lim^oo Sui (<?, to; t) where Su> m; i) is the time depen- 
dent tensorial density correlation function. The indices 
I, I' and m refer to expansion coefficients of a product of 
spherical harmonics Y^O) and Y™ (Q')- 

Since we have 

used the q-frame representation for the correlators, i.e. 
q has been chosen along the z-axis, the correlators are 
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diagonal in m and m! and are real || . The reader should 
also note that the existence of a rotational axis of sym- 
metry of the species of both systems allows to use the 
spherical harmonics, only. Using the results from ref. || 
and jn],[ll| one finds the molecular MCT-equations for 
these parameters: 

F(<7, m) = 

[S-\q } m) + m)K(q, m; {F( ? , m)})S- 1 (q, m)] _1 (1) 

where 

the matrix elements of the functional /C(q, m; F(g, m)) 
are given by: 

K(q,m;{F(q,m)}) = 

J2 E ?r(9)([mfem;{F( g ,m)})]-T"'^'(g) (2) 

a=T,Ra'=T,R 

with the MCT-polynomial: 

(m(q, m; {F(q, m)}))?" = m^? (q, m; {F(q, m)}) 
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and 
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g , « = T 
yW+T) , a = R 



(3) 



(4) 



The explicit expressions for the vertices V a - a for arbi- 
trary q can be found in ref. Q and for the q-frame in 
ref. |I(]JllJ • The vertices V aa only depend on the static 
correlators Su> (q, m) and the direct correlation functions 
cw (q, m) which are related to each other by the Ornstein- 
Zernike equation: 



S(q,m) 

where p is the number density. 
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III. COMPARISON 

For two particular systems the above equations have 
been solved. The first one is the diatomic linear A-B 
molecule of ref. H were the static correlators had been 
obtained by a molecular dynamics simulation. The pa- 
rameters of the Lennard- Jones potentials had been cho- 
sen such that the diameter of atom A and B were equal to 
1 and 0.95 (in Lennard- Jones units), respectively. Such 
an asymmetry is mainly needed to prevent crystalliza- 
tion. The distance between the two centers of the atoms 
was 0.5 . 

The second system is a fluid of prolate hard ellipsoids of 



revolution. Here the aspect ratio Xo = | which is the 
ratio between the major axis a and the minor axis b of 
the ellipsoids has to be chosen. If one wants to model the 
Lennard-Jones molecule by a hard ellipsoid the choice of 
the Lennard-Jones radii implies that Xo = 1.5. Note 
that Xo is not a fitting parameter. We also stress that 
the A-B molecules do not have a head-tail symmetry, in 
contrast to the ellipsoids. The value for Xo is in a pa- 
rameter range where no ordering due to a nematic phase 
is expected ]l2| , p^[ . The static correlators are obtained 
using Percus-Yevick (PY) theory. Using PY theory au- 
tomatically prevents crystallization since it is unable to 
yield a periodically ordered solution. PY, however, is in 
principle suited to study an orientational transition like 
an instability due to a nematic phase Jl^l . 
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FIG. 1. Static structure factor for the A-B Lennard-Jones 
molecule (solid lines) at the glass transition T c =0.475 ob- 
tained from ref. |^| . The components of the static structure 
factor for the hard ellipsoid of revolution with an aspect ra- 
tio of Xq — 1.5 and with, at the glass transition, a packing 
fraction of 0= 0.575 are plotted with dashed lines. From up 
to down the 1=L'=0, nr=0, the r=m=0, 1=2 and the 1=1'=2, 
m=0 components are shown. 

In fig. ^ we have plotted with solid lines the three 
static correlators with I = I' = 0; I = 0, 1' = 2 and 
I = I' = 2 for the Lennard-Jones system right at the 
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ideal glass transition (T c — 0.475 in L-J units and a den- 
sity of p c = 0.752) . The first one ( I = I' = 0) is the 
center of mass component. The second component with 
I = and V = 2 couples the two correlators for I = I' = 
and I = I' = 2, whereas the third one, the "quadrupolar" 
correlator gives information on the orientational order of 
the system. In fig. |]we have further plotted with dashed 
lines the same components for the system of hard ellip- 
soids at the critical density (packing fraction ip c = 0.575 
or density p c = -^f- — 0.732) of the glass transition. 
Note the remarkable similarity with the exception of the 
q — > behaviour where the simulation result has some 
shortcomings. 

These static correlators were used to solve the equations 
(0)"© which are truncated at l max — 2. As already men- 
tioned above, the A-B molecules do not possess head-tail 
symmetry. 
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FIG. 2. Non-ergodicity parameter as they arise from solv- 
ing the mode coupling equations in the limit of t — > oo for the 
A-B Lennard- Jones molecule (solid lines). The non-ergodicity 
parameters for the hard ellipsoids with an aspect ratio of 
Xo = 1.5 and, at the glass transition, a packing fraction <j>= 
0.575 are plotted with dashed lines. From up to down the 
1=1'=0, m=0, the P=m=0, 1=2 and the 1=1'=2, m=0 compo- 
nents are shown. 

Therefore the static correlators with I and/or V odd do 



not vanish, in contrast to those for the ellipsoids. Accord- 
ingly these correlators were also taken into account when 
solving the MCT-equations for the A-B molecules. The 
resulting critical non-ergodicity parameters Fu> (q, m) are 
plotted in fig. || for the A-B LJ molecule with solid lines 
for the ellipsoids with dashed lines. It is clearly visi- 
ble that the dominating contribution for the ergodicity 
breaking results from the first maximum (at q ~ 6.5 [— 1 ) 
of the structure factor of the center of mass component 
(/ = I' = 0). This is the manifestation of the cage effect. 

IV. CONCLUSION 

In conclusion, we have demonstrated that a system 
of hard ellipsoids of revolution is able to describe the 
static structure and the critical non-ergodicity param- 
eters (for I and I' even) of a system of super-cooled di- 
atomic molecules rather well, provided the aspect ratio is 
chosen properly. This result is not so obvious, since the 
static input into the MCT-equations for the diatomic 
molecules also includes the correlators with / and/or V 
odd. Those correlators either equal one or zero in case 
of the head-tail symmetric ellipsoids. For atoms with as- 
pect ratios Xq « 1.5 the cage effect determines the glass 
transition. This relation however will change close to a 
nematic instability. On one particular example we have 
shown that the detailed relation between orientational 
and center of mass components and the detailed shape 
and even the exact symmetry of the two-particle poten- 
tial seems only to have minor influence on the ergodicity 
breaking. Of course, it would be interesting to investigate 
how far more complex molecular systems can be modeled 
by hard ellipsoids. The full phase diagram for the glass 
transition of hard ellipsoids will be presented elsewhere 
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